Anisotropic non-Gaussianity from vector field perturbations 
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We suppose that a vector field perturbation causes part of the primordial curvature perturbation. 
The non-Gaussianity parameter /nl is then, in general, statistically anisotropic. We calculate its 
form and magnitude in the curvaton scenario and in the end-of-inflation scenario. We show that 
this anisotropy could easily be observable. 



Q>^ ■ I. INTRODUCTION 

o : 

] The primordial curvature perturbation ( provides one of the few windows available to the early Universe, and is the 
• subject of intense interest at present. It may be that C is practically Gaussian. Then its Fourier components will have 
J>~^' practically no correlation, except for the one required by the reality condition. The latter is defined by the two-point 
, ^ . correlator, specified by the spectrum 'P(^{k) where k is the wavenumber. On cosmological scales, the observed CMB 
^ ■ anisotropy gives an almost scale- independent value T'^ ~ (5 x 10~'^)^. 

^ i According to typical scenarios for the generation of (, the principle signal for non-Gaussianity would be the three 
point correlator, specified by the non-linearity parameter /nl- At present there is only an upper bound |/nl| ^ 100. 
I— — 1[ Over the next few years, the bound will go down to |/nl| < 5 or so if there is no detection. 

- During inflation, the vacuum fluctuation of each canonically normalized light scalar fleld becomes a classical per- 
turbation, with a nearly scale-independent and Gaussian spectrum. It is usually supposed that one or more of these 
Q , perturbations is responsible for the primordial curvature perturbation In this case, C is statistically homogeneous 
■ and isotropic, which means that its correlators are invariable under translations and rotations. Then 7^^(fc) depends 
c/2 . only on the magnitude of a momentum (wave vector) k, and /nlI^i, ^2, ^s) depends on the lengths of the sides of a 
, triangle. 

[ Most of the proposals for generating C from the scalar field perturbations belong to one of two broad classes, 
PsJ . which are distinguished by their prediction for /nl- If C is generated during single-field infiation one generally has 
^ i I/nl| ^ 10~^. More strongly, this bound applies to all single-field models in the squeezed configuration fci ~ fc2 » fcs. 
. If instead C is generated at or after the end of infiation, /nl becomes almost scale-independent. In this case the 
prediction for |/nl| is usually at least of order 1 and can be as big or bigger than the observational bound. 

Under the assumption of statistical isotropy, the constraints on /nl from current observation at the two-a level are 
[26] 



O 

oo 
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9 < fhT' < 111, - 151 < /nI"" < 253, (1) 



where the label 'local' can be taken to mean the squeezed configuration fci ~ fc2 ^ fcs and the label 'equil' can be 
taken to mean the equilateral configuration ki = k2 = k^. The first result might be regarded as weak evidence for a 
k>( , /nl**' ^ 1, which if confirmed would rule out the generation of C during single-field infiation. 
j_j ■ Recently, it has been suggested that the perturbation of some vector field may generate part or even all of the 
[ curvature perturbation. Such a thing is possible because the vacuum fiuctuation of a vector field can generate an 
almost scale-independent and Gaussian spectrum although the condition for that to happen is rather special (e.g. it 
doesn't happen for a canonically normalized light vector field.) The contribution to C from a vector field perturbation 
is statistically homogeneous but not in general statistically isotropic. Then the spectrum may depend on the direction 
of k, and the bispectrum may depend on the orientation of the triangle of vectors (fci, k2, fca). 

Observational bounds on to statistical anisotropy have not received much attention, and are not mentioned in the 
otherwise comprehensive analysis by the WMAP team [26]. As far as we know, the only available result [27] concerns 
the spectrum. It is parametrized in the form 

rd^)=\l+gid-kf]vrik), (2) 
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where k = |k| and the hats denote unit vectors. After taking account of all possible uncertainties, the conclusion from 
this study is g < 0.31 or so. We will adopt a bound g < 0.1 for definiteness. 

Using the SN formalism [1-4], one can write general formulas for the contribution of a vector field perturbation 
to the spectrum and the bispectrum. They were evaluated in Ref. [11] for the case that C is generated at the end of 
inflation, and in Ref. [5] for the case that C is generated by the curvaton mechanism. (Vector fleld inflation was also 
considered in Ref. [5] but we shall not consider it here.) It was found that the contribution to spectrum is of the 
form shown in Eq. (2). An analogous form for the contribution to /nl is to obtained in this work. These are the 
tree-level contributions. The 1-loop contribution to is given in Ref. [5] but the 1-loop contribution to /nl is not 
known at the time of writing. 

Prom the formula for , one sees that its statistical anisotropy could easily be as big or bigger than the observational 
bound. The purpose of this paper is to consider also /nl- We want to know if a vector field contribution to /nl could 
be big enough to observe, bearing in mind that its contribution to V(; should respect the observational bound. We 
answer this question in the affirmative, and give explicit formulas for the dependence of /nl on the orientation of the 
triangle. 



II. /nl including vector perturbations 



The evolution of the curvature perturbation ( on superhorizon scales is most readily described using a separate 
universe approach [1-4]. In a recent paper [5] the formalism was extended to take into account the possible statistical 
anisotropy in C, where it was shown that once one includes perturbations of the vector field, the resulting curvature 
perturbation can be calculated up to quadratic terms using the following equation 



C (x) = N^6cb + N^^SA, + ^N^^ {S^f + ^N^^SA, 



(3) 



where A'' is the number of e-folds of expansion of the unperturbed universe, the lower case roman letters denote spatial 
indices and Einstein summation over those indices is assumed. The derivatives of A^ with respect to the fields are 
denoted as 



dN 



N 



dN 



d(t>dAi 



and AT^ 



dA^dA, 



(4) 



To obtain Eq.(3) it was assumed that the anisotropy of the expansion of the Universe is negligible. Even for the vector 
fields having a considerable contribution to the curvature perturbation the isotropic expansion of the Universe can be 
achieved at least in several ways, for example with the oscillating massive vector field [6-8], using a triad of orthogonal 
vectors [9], a large number of identical randomly oriented vector fields [10] or if the contribution of the vector field(s) 
to the total energy density is negligible [11]. However, in general perturbations generated by vector fields will induce 
statistical anisotropy in Q. In Eq.(3) we assume only a single scalar field and a single vector field A^. 
We define the power spectrum and the bispectrum through the Fourier modes of C as: 



(C (ki) , C (k2)) ^ {2^f <5(3) (ki + k2) ^T'c (ki) , 



(5) 



(C (ki) C (k2) C (ks)) = (2^)' S^""^ (ki + k2 + kg) (kl, k2, kg) , 

where normalization of Fourier components is chosen to be such that 

C(k)^ / C(x)c-*'^-dx. 



(6) 



(7) 

Note that the power spectrum and the bispectrum are dependent on the direction of k. The bispectrum (ki , k2, kg) 
can be further separated into three parts: one due to perturbations in the scalar field, another part due to the vector 
field and a mixed term: 



(ki,k2,k3 



TJT:3V^{ki)V4,{k2) + C.p. 

-TgTgT'^ {ki)Mt (k2) + 5 perm. 



B^^(ki,k2,k3) = --A^^A^, 
Att'^ 

BA(ki,k2,k3) = -3-3X,(ki)A^,^A4,(k2)+C.p. 
K1K2 



(8) 
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where "c.p" stands for "cyclic permutations" and ki, k2, are the moduH of the vectors ki, k2 and ks. 

The power spectrum (fc) in the above equations depends only on the modulus of k because we assumed that the 
expansion during inflation is isotropic. The vector M i (k) characterizes perturbations of the vector field: 

M{\<i) = V+{k)NA N^+p(fc)k(k-N'^) +ig(fc)kx . (9) 

In this expression Na is the modulus of the vector N"^, N"^ and k are unit vectors defined by = 'N'^ /Na and 
k = k/A:. The power spectrum for the longitudinal component is denoted by Vq (k) while 7^+ (fc) and V- {k) are the 
parity conserving and violating power spectra defined by 

v± = ^{rR±rL), (10) 

with VR{k) and Vhik) denoting the power spectra for the transverse components with right-handed and left-handed 
polarizations. Also we have defined p{k) and q{k) as 

Vo - V+ V- 
P= — and q=——. (11) 

Because of the isotropic expansion during infiation the power spectra Vo and V± depend only on the modulus of k. 
The curvature perturbation power spectrum 'P^ (k) may be separated into isotropic and anisotropic parts [12]: 



Vc{k) = r'^°{k) 



1 



'(fc) (n^ -k)' 



(12) 



where the amount of anisotropy at each scale is parametrized by 



The isotropic part of the spectrum is 



g^Ni'^^. (13) 

V'^° = NlV^ + NlP+. (14) 

In this paper we will calculate /nl for both configurations: equilateral, in which ki = k2 = k^, and squeezed, in 
which fci ~ fc2 » fca- In the equilateral configuration the bispectra from Eqs. (8) become 

e;'i""(ki,k2,k3) - 3NlN^^Pl{k,), 

S^7'(ki,k2,k3) = ^N^N^,^V4,{k,)[M^(k,)+M,{k2)+M,{k3)], (15) 

e^'l""(ki,k2,k3) = X,(ki)iV,^X,(k2)+C.p., 

where we have defined for the equilateral configuration 

S^°^""(ki,k2,k3)^ (^|^yi3°'^""(ki,k2,k3), (16) 

and g^'i"'' — goquii _|_ _|_ In this case the non-linearity parameter /^l"'' is expressed through the power 

spectrum and the bispectrum as: 

5^NL 3piso(fc)2 ■ 

Observations give a limit on the anisotropy g ^ 0.1 [27]. Therefore, since the anisotropic contribution to the curvature 
perturbation is subdominant compared to the isotropic one, we have included only 7',^''° into the above expression 

Of/Nf'- 

For the squeezed configuration two of the vectors have almost identical lengths but opposite directions, ki ~ — k2, 
but the third vector k3 is of much smaller modulus than the other two and almost perpendicular to them. For this 
configuration Eqs. (8) take the form 

S;°^'^'(ki,k2,k3) = 2N'^N^^V^{ki)r^{k3), 

B^°5^'(ki,k2,k3) = -iV^iV^,;{7'^(fci)A^,(k3)+7'^(fc3)Rc[X,(ki)]}, (18) 
S^^^'(ki,k2,k3) = 2Rc[A<,(ki)] A^4Rc[A^, (k3)], 



where Re [. . .] means the real part and M"'^'*' (ki, k2, ks) is defined similarly to Eq. (16) 



(ki,k2,k3) 




ilocal 



(ki,k2,k3) . 



(19) 



Then, the nonlinearity parameter /j 



^*local 
NL 



in the squeezed configuration becomes 




plocal 




(20) 



Having defined our notation, in the following sections we will consider two particular examples. But before that, 
we wish to point out an important subtlety concerning the vector field A^. By Ai we refer to the spatial components 
of the physical vector field Ai = Bi/a, where a is the scale factor in FRW universe. The field Bi, which enters the 
Lagrangian, is the comoving vector field with the expansion of the universe factored out. In FRW spacetime the 
temporal components of the physical and comoving fields are the same, = Bq [6, 7]. 



In this section we study the non-Gaussianity in the curvature perturbation for a model of slow roll infiation 
with an additional U (1) vector field which decays some time after reheating and contributes to the total curvature 
perturbation following the curvaton mechanism [14]. For this contribution to be non-negligible, the vector field must 
undergo particle production during infiation and obtain a superhorizon spectrum of perturbations. But the massless 
U (1) field is conformally invariant and, consequently, its quantum fiuctuations are not amplified during infiation. 
This means that in order for such a vector field to undergo particle production we have to brake its conformality [15]. 
Another problem is that a dominant homogeneous (homogenized by infiation) vector field could make the expansion 
of the Universe strongly anisotropic, which is in contradiction with observations. 

One way out of this problem can be the curvaton mechanism [14, 16, 17]. The usual curvaton scenario incorporates 
two scalar fields: one that drives infiation and another one, called curvaton, which produces the curvature perturbation. 
During infiation the curvaton is subdominant. However, after reheating the universe is radiation dominated and its 
energy density is diluted as pr oc a~^. If the energy density of the curvaton field decreases slower than a~^, at some 
moment it can dominate (or nearly dominate) the Universe and impose its own curvature perturbation. This is the 
basic idea of the curvaton mechanism. 

Here we consider a massive vector field acting as the curvaton. Before dominating the vector curvaton field is rapidly 
oscillating in a quasi-harmonic manner. As shown in Ref. [6], the oscillating vector field behaves as a pressureless 
isotropic fiuid and can dominate without generating a large-scale anisotropy. 



Here we obtain analytic expressions for the non-linearity parameter /nl without assuming a specific vector curvaton 
model. In contrast to the original curvaton idea we include as well perturbations generated during infiation by the 
light scalar field (for a similar study in the scalar curvaton case see Refs. [18, 19]). 

When some time after reheating the mass of the vector field becomes bigger than the Hubble parameter, the field 
starts to oscillate. In Ref. [5] it was shown that for oscillating vector field and are equal to 



III. ANISOTROPY IN THE VECTOR CURVATON MODEL 



A. The generic treatment 



2 A 

s^'a^' 



(21) 




(22) 



where A = |A| is evaluated just before the vector field decays and the parameter r is defined as 
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with PA being the energy density of the vector field just before its decay (taken to be sudden), Q,a = Pa/ P and 
p = PA+ Pr- Using Eq. (21) the isotropic part of the total power spectrum in Eq.(14) becomes 



Vf = NlV^[\-.(i'^\, (24) 



where we defined 



Then the vector part of the bispectrum for equilateral configuration in Eq.(15) reduces to 

BT'\]^^MM) = (Yx) ^n(A:i)P+(fc2)(l+p(fcl)A2+p(fc2)A2+AiA2 



q{ki)q{k2) - ^p{ki)p{k2) 



+t^^-iAl+ A1A2 + Al) [Aip (fci) g (fca) - A^p (fca) q (ki)] + (fci) q (fcs) | + c.p. (26) 

In the above we used the notation ^1 = A • ki etc., where A = A/ A. Because the configuration of wavevectors ki, 
k2 and ka is equilateral, with the angle between any two of them being 27r/3, we find ki • k2 = ki • ka = k2 • ka = — ^■ 
Eq.(26) simplifies further if we consider a scale invariant power spectrum and the expression for /^l"'' becomes: 



5/NL -P-P^^, iv, + pv^f ' ^ 



± 



where we have taken into account that the non-Gaussianity generated during the single field infiation is negligible. 
The quantity A± < 1 is the modulus of the projection of the unit vector A onto the plane containing the three vectors 
ki, k2 and k^. The calculation of A± in the equilateral configuration is explained in more detail in the Appendix. 
For the squeezed configuration the bispectrum from the vector field perturbation in Eqs. (18) becomes 



2 r \^ 1 



S;?-'(ki,k2,k3) - 2(^-j jV+{k,)V+{k3)[l+p{k,)Al+p{k:,)Al]. (28) 



Working as in the equilateral case, we find that the non-linearity parameter for the scale invariant power spectra is 

6 .local _ o2p2 _^ l+pA\ 



In this equation ip is the angle between the vectors ki and Aj_. 

As one can see from the above equations, /nl is, in general, dependent on A±^, in both configurations. This means 
that /nl is anisotropic and that the amount of non-Gaussianity is correlated with the statistical anisotropy. However, 
from Eqs. (10), (11) and (13) it is clear that, if particle production is isotropic (i.e. Vq = 7-"+ and V- = 0) then 
p = q = Q and the above expressions for /^l"'' ^-nd /nl**' become isotropic too and both reduce to /nl = 5/4r as 
in the scalar curvaton scenario, where we have assumed that ^ 7^+, i.e. that the dominant contribution to the 
curvature perturbation is due to the vector curvaton field only. 

B. /nl for non-minimally coupled vector curvaton 

In Ref. [6] it was shown that a vector field can attain a scale invariant perturbation spectrum if it's mass during 
infiation is equal to m? = —2H^. One way to achieve a negative mass squared of this magnitude is to introduce a 
non-minimal coupling of the vector field to gravity of the form ^RB^B^, where R is the Ricci scalar. The idea of 
such a non-minimally coupled vector curvaton was introduced in Ref. [8]. In that paper and in Ref. [5] it was shown 
that the power spectra for different polarizations are 



n=(^ , and 7^0 = 2 -) . (30) 
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One notices that the parity conserving transverse power spectrum and the power spectrum generated during the 
single scalar field infiation are equal, i.e. V+ = V^. Thus the isotropic part of the curvature perturbation spectrum 
can be written as 

Vf°^V^Nl{l + p). (31) 
While the anisotropy parameter from Eq.(13) becomes 

.-TT^. (32) 

Using Eq.(30) we find 

p=\ and g = 0. (33) 



Thus, the anisotropy in the vector field is rather strong, which means that it will have to remain subdominant, i.e. 

j'Cqui 



< 1. Using this and Eq. (27), the f^^'^ for the non-minimally coupled vector curvaton is found to be 



Similarly, /^l^' for the squeezed configuration in Eq. (29) is 



2 



^/;r = 2|-(1 + Ai) (35) 



Since V+ = \Vo =1^^ = (|^) , for the typical values of the perturbations we have 5(j) ^ SAi ~ H. This means that, 
in order for the vector field contribution to be subdominant, we require Na ^ (c.f. Eq. (3)). Hence, /3 <C 1 and 
g — 13. Thus, in view of Eqs. (34) and (35), we see that /nl ^ 9^ I^A- Therefore, we find that the non-Gaussianity is 
determined by the magnitude of the statistical anisotropy. 

This prediction is valid in the regime <C 1 which corresponds to ^\ > V^P, which impHes /nl ^ 9^^^/ 

For smaller fi^, the contribution of the vector field perturbation to C is of order 17^[i5A/(i5A2)^/^]. In other words, it 
is of order and is the square of a Gaussian quantity. The resulting prediction for its contribution to /nl would be 
given by a 1-loop formula which has not been evaluated at the time of writing. 



IV. ANISOTROPY GENERATED AT THE END OF INFLATION 



As another example, we consider the generation of an anisotropic power spectrum at the end of infiation. The idea 
is based on Ref. [21] where it was shown that in hybrid infiation models the generation of curvature perturbations can 
be realized due to inhomogeneous end of infiation. Yokoyama and Soda [11] used this idea to generate the anisotropic 
contribution to the total curvature perturbations. In their model the anisotropy is generated at the end of infiation 
due to the vector field coupHng with the waterfall field. In this section we calculate the non-Gaussianity of the model 
in Ref. [11] using our notation in section II. 

In this model there are two components of the curvature perturbation: one generated during infiation and an 
anisotropic one, generated by a vector field at the end of infiation: 

C = Cinf + Ccnd- (36) 

The first component is due to the perturbation generated during infiation. The second component is due to the 
perturbation of the vector field. Ref. [11] considers a massless, U (1) vector field. Without parity violating terms the 
power spectra for left handed and right handed polarizations are equal, while the longitudinal polarization is absent 
for a massless field. In this situation we find that parameters p (fc) and q (fc) from Eq.(ll) become 

p = -l and 9 = 0. (37) 

The conformal invariance of the (7(1) vector field is broken through a non-canonical kinetic function of the form 
/^ (t) Ffj^i.F'^'^ , where = d/^iBi, — d^B is the field strength tensor and - comoving vector field. This form of 
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conformal invariance braking was considered in many papers (e.g. Refs. [7, 22-25]) where it was found that a scale 
invariant perturbation spectrum is obtained if / oc a: 



So the isotropic part in Eq.(14) of the power spectrum becomes 

This is of the same form as with the vector curvaton model but with different (3: 

Na ■ ' 



The anisotropy parameter in Eq.(13) in this model becomes 



f3 



1 + 13 



Taking into account Eq.(37), the vector Mi (k) in Eq.(9) reduces to the simple form 



k • k 



M (k) = NaV^/-^ [N' 
To calculate /nl we consider a specific example of hybrid inflation with the potential 
V 



X, S^) = Vo + ^m202 _ 1 ^2^2 ^ t.Xx^ + ^X^^^x^ + ^XAx'Bf^B, 



where 4> is the inflaton and x is the waterfall field. The effective mass of the waterfall field for this potential is 



'off 



-m^ + X^cfP' - XAAiAi, 



(38) 



(39) 



(40) 



(41) 



(42) 



(43) 



(44) 



where Einstein summation is assumed and we used Ai = Bi/a and the Coulomb gauge in which Aq =Q and diA^ = 0. 
Infiation ends when the infiaton reaches a critical value 4>c where the effective mass of the waterfall field becomes 
tachyonic. But one can see from Eq.(44) that the critical value is a function of the vector field 0c ~ 4'c (A). With this 
in mind the vectors and can be readily calculated: 



and 



where we have defined 



ON 92 



dN d(j)c 
dSr dAi 



XaA 

X^ <j>c ' 



DAW A3 



N, 



OA, dAj 4>cN, 



^ and C^,/^^, 
d(j)c y Xa A' 



AiAj 1 , 



(45) 



(46) 



(47) 



where A is evaluated at the end of infiation and we used the fact that Ncc/N^ ^ N^^/N"^ ^ O (e) under the slow 

roll approximation [21], where e is the slow roll parameter defined as e = \M'^ {V /V)^ , with the prime denoting 
derivatives with respect to the infiaton. As mentioned earlier the total of perturbations consists of two components: 
perturbations of the scalar and vector fields. This gives the following bispectrum in the equilateral configuration 



S^'^""(ki,k2,k3) = 67"(ki,k2,k3)+e^^"''(ki,k2,k3) = 

= WlNjN^^ + [M, (ki) N^^M, (k2) + c.p.] = 



'Nc(l)c 



iC'-l)-[-,C^-l]Al-^^A 
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The mixed term S^X'' absent from Eq.(48) because in this model iV^ = 0. By inserting Eq.(39) into (48) we obtain 



6 j-cquil 2 



{c- - 1) - -c 



1 Ui 



16 



(49) 



where the slow parameter 77 is equal to = m^Mp/Vo and MpNc = l/^/2e^, with Cc being the slow roll parameter at 
the end of inflation. Similarly, for the squeezed configuration we find 



6 j-ic 



{C'-l){l-Al)-'-{sm^)'Ai 



(50) 



Again, we find that /^l"'' ^ind /nl**' are functions of A±, i.e. they are anisotropic and correlated with the statistical 
anisotropy. Also the level of non-Gaussianity is proportional to the anisotropy parameter squared, /nl oc g^, as in 
the non-minimally coupled vector curvaton model. 



V. DISCUSSION 



In this paper we have calculated the amount of the non-Gaussianity generated by the anisotropic part of the 
curvature perturbation. We have also considered two specific models to generate the anisotropic curvature perturbation 
from the vector fields and calculated the non-Gaussianity in detail for those models. The results were given for the 
equilateral and squeezed configurations in Eqs. (34) and (35) for the model of non-minimally coupled vector curvaton 
and in Eqs. (49), (50) for the end-of-infiation scenario. We have shown that /nl generated by the vector field is 
anisotropic and that it is correlated with the amount and direction of the statistical anisotropy. 

Although in some specific models it would be possible to generate the statistically isotropic curvature perturbation 
only from the vector field, however, in general, perturbations generated by the vector field are anisotropic. From 
observations the bound on statistically anisotropic contribution to the total power spectrum is constrained to be less 
than about 10 %. In this case one can estimate the maximal /nl caused by the statistically anisotropic curvature 
perturbation on fairly general grounds. 

Let us assume that the non-Gaussianity is produced solely due to the vector field perturbations. If this is so 
then Eqs. (17) and (20) suggest /nl ^ Bc^/Vq, where we consider that the anisotropic contribution to the curvature 
perturbation is subdominant, i.e. Vq ~ 7'^''°. According to Eqs. (16) and (19) we have ~ M'^Naa, where M = | Al| 

and Naa = 11^^11- Thus, we have /nl ^ M^Naa/V^ 

Now, M. depends on the mechanism which breaks the conformality of the vector field and is responsible for the 
generation of its superhorizon perturbation spectrum. If this mechanism does not introduce additional mass scales, 
then, on dimensional grounds, we expect the anisotropy in the vector field perturbation to be of order unity, i.e. 
\p\, \q\ ^ 0(1), barring cancellations such as due to parity invariance (which results in g = 0) or an isotropic particle 
production, which gives p = q — and generates no statistical anisotropy. In our examples in Sees. Ill and IV we indeed 
analysed such a situation, where the effective mass-squared of the vector field during inflation = = — 2i7^ or 

the time dependence of the kinetic function f / f = H were both determined by the dynamics of the expansion and 
given by H, the only scale in the theory. If the anisotropy in the vector fleld perturbation is of order unity, then 
Eq. (9) gives M ^ VaNa, where Va = 27^+ + is the power spectrum of the total vector fleld perturbation given 
by Va = X^i I'^^iP in the superhorizon limit. Putting the above together we obtain 

JNL ~ ^ ~ p2 ~ ^2 ■ y^^) 

Since we are working in the regime where <C 1 we expect the higher order contribution to C in Eq. (3) 

from the vector fleld to be subdominant, i.e. Na^A > Naa^A^. Considering that the typical value of the vector fleld 
perturbation is 5A ^ VVa, we obtain the bound /nl < {V]^^ NaY /V^. As evident from Eq. (3), the contribution of 

the vector fleld to C is given by Ca ^ \/T\a ^ Na^/Va, where T^^a i^ the power spectrum of the anisotropic curvature 
perturbation. With this in mind, the upper bound to /nl becomes 

hi^T'f<{^)"' ■ (52) 

Because the vector fleld contribution to the total curvature perturbation must be subdominant, Eq. (12) suggests that 
the anisotropy of the curvature perturbation is 5 ^ ^Ca/^C- Using this and also that Vq ~ b'K 10~^, we flnd that the 
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maximum value of /nl generated by the statistically anisotropic contribution to the curvature perturbation has to be 



It is clear that the study of vector field contributions to the primordial curvature perturbation is just beginning. 
Given a scale invariant and Gaussian vector field perturbation, the calculation of C from the SN formalism is straight- 
forward and should now be done for the full range of scenarios that have already been explored for the contribution of 
scalar field perturbations. Also, the Feynman graph formalism available for the scalar field case should be generalized 
to cover the vector field contributions. At a deeper level, one also wishes to understand how the scale invariant 
perturbation can be generated. One would also like a fuller understanding of the generation of perturbations from 
the vacuum fiuctuation when the expansion of the unperturbed Universe is anisotropic, since that can easily happen 
in the presence of vector fields. 

For the moment though, it is most urgent to confront specific predictions for the form of the anisotropy with 
observation. Specifically, we want to know what constraint is placed by observation on a contribution to /nl of the 
form in Eqs. (34), (35) or Eqs. (49), (50). It seems quite possible that with the latter prediction, valid for the squeezed 
configuration, one might find a nonzero value at better than the 2-(t level that is already found for the isotropic case. 
Such a finding, if confirmed, would be a smoking gun for a vector field contribution to the curvature perturbation. 
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Our results in Eqs. (34), (35), (49) and (50) apply if /nl is below this value. 
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APPENDIX 

CALCULATION OF A± IN THE EQUILATERAL CONFIGURATION 

First note that in the equilateral configuration ki + k2 = — ka. This gives Ai + A2 ~ — A3 and 

Al + Al + Al = 2 {Al + A1A2 + Al) ■ 

A1A2 + + = -{Al+ A1A2 + Al) ; (54) 
AlAl + AlAl + AlAl = {Al + A^A2 + Alf . 

Let us define a vector A_l which is the projection of A to the plane containing vectors ki, k2 and ks. Then the scalar 
product of these vectors and A is the same as the product with Aj^: 

A • ka = • ka, (55) 

where a = 1, 2, 3. 

Without loss of generality we can assume that the angle between Aj^ and ki is ip: 

= A • ki = A^ • ki = A±^ cos (p, (56) 

where A±_ = |A_l|. In equilateral configuration the angle between vectors ki and k2 is 27r/3, and A2 becomes 

< / 27r\ , /l V3 \ 

A2 = Aj^ • k2 = A^ cos I + — I = - Ax I - cos + — sm 1 . (57) 



From the last two equations we get 



Al+AiA2+Al = -^A\. (58) 



Putting this result back into Eq.(54) we find 

Al + Al + Al = lAl- 



A1A2+A2A3 + A3A1 = -|Ai; (59) 

42 /1 2 I A2 a2 , 4242 _ _9_44 
12 ^ 2 3 ^ 3 1 — T6 -L ' 



